Pulsar Kicks With Sterile Neutrinos and Landau Levels 
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We use a model with two sterile neutrinos obtained by fits to the MiniBoone and LSND ex- 
periments. Using formulations with neutrinos created by URCA processes in a strong magnetic 
field, so the lowest Landau level has a sizable probability, we find that with known paramenters the 
asymmetric sterile neutrino emissivity might account for large pulsar kicks. 
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I. INTRODUCTION 

The gravitational collapse of a massive star often 
leads to the formation of a neutron star, a pulsar. 
It has been observed that many pulsars move with 
linear velocities of 1000 km/s or greater. See Ref.[l[ 
for a review. We have investigated the pulsar kicks 
which arise from the modified URCA processes in 
the time interval 10-20 sec after the supernova col- 
lapses, with a strong magnetic field and temperature 
so that the population of the lowest Landau level is 
approximately 0.4 of the total occupation, and we 
find large pulsar kicks 0. 

The largest neutrino emission after the supernova 
collapse takes place during the first 10 seconds, with 
URCA processes dominant. The possibility of pul- 
sar kicks from anisotropic neutrino emission due to 
strong magnetic fields during this time was discussed 
more than two decades agoQ. It has been shown Q 
that, with the strength of the magnetic field ex- 
pected during this period, the lowest Landau level 
has a sizable occupation probability, which produces 
the neutrino emission asymmetry that is needed for 
pulsar kicks. However, due to the high opacity for 
standard model neutrinos in the dense region within 
the neutrinosphere, few neutrinos arc emitted, and 
the pulsar kick is not obtained^. Sterile neutrinos 
with a small mixing angle have small opacities. It 
has been shown @ that using the model of Ref.[l] 
and assuming the existence of a heavy sterile neu- 
trino (mass > 1 kev), with mass and mixing angle 
constrained to fit dark matter, the pulsar kicks might 
be explained. 

Recently, the MiniBooNE Collaboration found 
that the data for electron neutrino appearance 



showed an excess at low energies, in comparison to 
what was expected in the standard modelfji]. This 
data, along with the LSND data, has been analyzed 
in a model with two light sterile neutrinos Q, and 
compared to MiniBooNE data@. The mixing angles 
of two light sterile neutrinos were extracted. See, 
however, Ref . [ic| . which questions the accuracy of 
the results of Ref [1]. 

In the present paper we use the fits of Refs.Jl, Q 
with two sterile neutrinos to investigate the possi- 
bility of obtaining the large pulsar velocities which 
have been observed. As we shall show, our model 
differs from that of Ref[6j in that with a much larger 
mixing angle there is a higher probability of ster- 
ile neutrinos, but a much smaller effective volume, 
due to a larger opacity. However, as we shall show, 
since the mean free path is much larger than those 
of standard neutrinos, under the conditions in which 
standard neutrinos produce a pulsar velocity of 2- 
300 km/s, the MiniBoone/LSND sterile neutrinos 
can give a kick of more than 1000 km/s. 



II. ASYMMETRIC STERILE NEUTRINO 
EMISSIVITY AND PULSAR KICKS IN 
LIGHT TWO-STERILE NEUTRINO MODEL 

Within about 1 second after the gravitational col- 
lapse of a large star, the neutrinosphere with a radius 
of about 40 km, with temperature equilibrium, is 
formed. For about 10 seconds about 98% of neutrino 
emission occurs, with neutrinos produced mainly by 
URCA processes. Due to the strong magnetic field, 
neutrino momentum asymmetry is produced within 
the neutrinosphere, but with a small mean free path 
they are emitted only from a small surface layer of 
the neutrinosphere, and the pulsar kick cannot be 
accounted for. If a standard active neutrino, say the 
electron neutrino, oscillates into a sterile neutrino, 
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it will escape from the protoneutrino star and neu- 
trinosphere, unless it oscillates back into the active 
neutrino. The mixing angle plays a key role. In the 
work of Fuller et alp] the mixing angle is so small 
that the sterile neutrinos are emitted. In the present 
work the starting point is the analysis of MiniBooNE 
and LSND data, with the two or more sterile neu- 
trinos with small masses and large mixing angles. 
Before we can proceed, however, it is essential to 
determine possible effects of the high density and 
temperature of the medium on the mixing angles. 

A. Mixing Angle in Neutrinosphere Matter 

It has long been known that dense matter can ef- 
fect neutrino states. The famous MSW effect (IllIT^ 
for understanding solar neutrinos, and the study of 
oscillations of high energy neutrinos [l3j are studies 
of mixing of active neutrinos in matter. There have 
been many other studies. In the present work we are 
dealing with sterile/active neutrino mixing given by 
the mixing angle 9 m in neutrinosphere matter 



a <1.0, p=b MeV, and {5m) 2 = 1.0 ev 2 @, dc}, we 
find 

2nV T 

~ 5.1 x l(T 3 & 2 a 4 < cos(29) . (4) 

(dm)- 1 

Therefore the mixing angle in the neutrinosphere 
medium is approximately the same the vacuum mix- 
ing angle. This agrees with Ref@. 



B. Emissivity With a Light Sterile Neutrino 

We now use the fits to MiniBooNe and LSND with 
light sterile neutrinos to estimate pulsar kicks. The 
MiniBooNE results are consistent with the LSND re- 
sults only if there are at least two sterile neutrinos. 
Models with three sterile neutrinos have also been 
considered [Tol liH]. Fits to the MiniBooNE experi- 
ment and the LSND results by Ref.[8] in Ref.[|[ with 
two sterile neutrinos are shown in Fig. 1. 



\vi > = cos9 m \v e > -sin6 m \v s > (1) 
\i> 2 > = sin9 m \v e > +cos9 m \v s > . 

In the work Ref[6] it was shown that the mixing 
angle for sterile neutrinos that can account for dark 
matter as well as those produced in the neutron star 
core is almost the same as the vacuum value. Start- 
ing from the much larger mixing angles for the sterile 
neutrinos that seem to account for the MiniBoonE, 
LSND data, we need the value of the mixing angles 
in the neutrinosphere, as we discuss below. The ef- 
fective mixing angle in matter, 9 m can be related to 
the vacuum mixing angle, 9 bvfhlj 



sin 2 (29 rl 



sin 2 (29) 



sin 2 (29) + (cos{29) 



2pV T ^ 2 
(8m)*) 



•(2) 



In Eq@ V T is the finite temperature potential, 
while the finite density potential due to asym- 
metries in weakly interacting particles has been 
dropped as it vanishes when temperature equilib- 
rium is reached[14]. A convenient form for V T , with 
the background of both neutrinos and electrons in- 
cluded, is given in Ref[l|| 



45a 



2 9 w (l + 0-5cos 2 9 w ) P T 4 , (3) 



with Gf,9w the standard weak interaction param- 
eters and a = 1/137. Assuming T=20a MeV, with 
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FIG. 1: sm 2 (26»i s )=0.004; sin 2 (26» 2s )=0.2 
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From the Sterbenz/Maltoni-Schwetz fits one finds 
for the mixing angles of the two sterile neutrinos: 



(sin29 ls ) 2 
(sin28 2s ) 2 



0.004 
0.2 , 



(5) 



and the masses are negligibly small. Note that this 
is in contrast to the parameters of Ref . Q , with the 
constraint of dark matter giving a mixing angle of 
(sin26dm) 2 — 10 -8 , and a mass greater than 1 keV. 
Note that in a recent analysis by Maltoni (TFJ] it 
was found that the fit [8| to the MiniBooNE and 
LSND data had problems with disappearance data, 
although there are rather large error bars. In our 
present work we will use values for (sin26) 2 in the 
range 0.2 to 0.004 to estimate the pulsar kick. 

The probability of asymmetric emission, giving a 
pulsar kick, does not depend directly on the sterile 
neutrino mass in our model, but is proportional to 
the (sin26 s ) 2 . It is the large mixing angles found in 
fits to MiniBoone and LSND that lead us to carry 
out the investigation in the present paper. 



C. General Formulation of Neutrino 
Emissivity With a Strong Magetic Field 

The neutrino emissivity is give n in general form 
in many papers, e.g., see Refs |l7t IT^|: 



n 
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d 3 q e 



(2tt) 3 2c^(27r) 3 J (2tt) 3 
(2tt) 4 £ J_ w ^MtM (6) 
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where M is the matrix element for the URCA pro- 
cess and T is the product of the initial and final 
Fermi-Dirac functions corresponding to the temper- 
ature and density of the medium. The main source 
of the asymetric emissivity that produces the pul- 
sar velocity is the fact that the electron has a large 
probability to be in the lowest (n=0) Landau level. 
See Refs [19|, [2(3] for a discussion of Landau levels. 
The asymetric emissivity can be seen by considering 
the weak axial intereaction, Wa, 



W A 



G 



9Axt l ■ <?Xr 



(7) 



^(q e ), ^(q l/ ) ) where q e and q" are the electron and 
antincutrino momenta, respectively. The key to the 
asymmetric emission is given be the trace over the 
leptonic currents, Tr[ljlj], 



d 2 q<\Tr[l\l ] ~ 8KE e [(qy5 i3 + (<f)^ 



J3 



-<W) 3 K<f = B = z) 



(8) 



with the magnetic field B in the z direction. We only 
consider the weak axial force, which is dominant. 
Using the relationship given in Eq® one can show 
that the result of the traces and integrals over the 
axial product matrix element has the form (B = z) 



\m a \ 2 cx {<ry 



(9) 



Details are given in Ref [2J, where it is shown that 
the asymmetric neutrino emissivity, using the gen- 
eral formulation of Ref [171. is 



0.64 x 10^T 9 7 P(0) x / 



erg cm 3 s 1 = p ns c 



Veff^t 



(10) 



where Tg = T/(10 9 K), p ns is the neutron star mo- 
mentum, P(0) is the probability of the electron pro- 
duced with the antineutrino being in the lowest Lan- 
dau state, f=.52 is the probability of the neutrino 
being at the + z neutrinosphere surface 0, V e ff 
is the volume at the surface of the neutrinosphere 
from which neutrinos are emitted, and At ~ 10s is 
the time interval for the emission. 

We derive P(0) and V e //, the effective volume for 
the emissivity, in the next two subsections. 



D. P(0) = Probability for the Electron to be in 
the n=0 Landau Level 

Just as in our previous work in which Landau 
levels play a crucial role[2|, only the lowest Lan- 
dau level, for which the helicity is -1/2 (rather than 
±1/2 as with the usual Dirac spinors) gives asym- 
metric emission. The probability that an electron in 
a strong magnetic field is in the lowest (n=0) Landau 
level, P(0), can be calculated from the temperature, 
T, and the energy spectrum of Landau levels [l9l.[20T|. 
A particle with momentum p and effective mass m* 
in a magnetic field B in the nth Landau level has the 
energy 



with G 



fJA 



1.26, the x are the nu- 



cleon spinors, and the lepton wave functions are 



E L (p,n) 



o 2 + (m*) 2 + 2(m* e ) 2 ^-n , (11) 
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with B c = 4 x 10 13 G, and m* is the effective mass of 
the electron at the high density of the protoneutron 
star and neutrinosphere. 

From standard thermodynamics the probability of 
occupation of the n=0 Landau level, P(0), is given 
by0: 



x(l + e (^-^)/ feT ) 



(14) 



P(0) 



F(0) 



(12) 



where F(n), with magnetic field B, temperature T, 
and chemical potential /i, is 



F(n) 



dp 



m p -E L (p, n)} 2 



l + exp[(E L (p,n)- n)/T] 



(13) 



The electron energy is restricted to magnitudes 
greater than /j,, but the integrals in Eq.([7]) are in- 
sensitive to Pmin, so we take p m in = as in Ref@. 

We agree with the estimate of Ref.[(| for P(0). 
Note that if we had used the free electron mass, m e , 
in the Landau energies (Eq.©) we would have ob- 
tained a much smaller value for P(0). For B= 10 16 
G, fi = 40 MeV, m*=4 MeV and T^ sphere = 20 
MeV, P(0) ~ 0.3. This is similar to our estimate 
of P(n=0) ~ 0.4 at the surface of the protoneutron 
star at about 10 seconds[2j. Therefore our result for 
asymmetric emissivity differs from that of Fuller et 
al [|| mainly in that we have a much larger mixing 
angle, and a much smaller effective volume, since the 
sterile neutrinos oscillate back to active neutrinos 
within the neutrinosphere; and therefore our emis- 
sion only takes place near the surface of the neu- 
trinoshere. However, in contrast with purely active 
neutrino emission in which the opacity results in very 
small pulsar kicks [H, the sterile neutrinos have a 
much larger effective volume, and can therefore pro- 
duce much larger pulsar velocities 



E. Estimate of V e ff= Effective Volume for 
Emission 

To estimate V e ff we make use of the early study 
of opacity in about the first 20s of the creation of 
a neutron star via a supernova collapse [U [23|, [24[ , 
and a recent detailed study of neutrino mean free 
paths[22j. Since the mean free path of the sterile 
neutrino is determined by that of the standard neu- 
trino to which it oscillates, A, we make use of studies 
of active neutrino mean free paths. First note that 
the neutrino mean free path is given by 



1/A = 



d 3 p 
(2tt) 3 



M fi [l-n(q)] 



where Mfi is the weak matrix element and n(q) is 
the Fermi distribution. For the calculation of the 
sterile neutrino, for which Mfi — 0, one can use 
the value of 1/A with a factor of sin 2 (29) from the 
matrix element and another such factor from the 
occupation probability. From the results of previous 
authors, for T in the 10 to 20 Mev range and /i in 
the 20 to 40MeV range, we estimate that A ~ 1.0cm 
This gives a range for the effective sterile neutrino 
mean free path 



A, ~ 5.0 to 250 cm 



(15) 



For a neutrinosphere radius of 40 km, with A s < < 
R v this gives us V eff = (47r/3)(i? 3 - (R v - A s ) 3 ) ~ 

From Ea.lflD). R u = 40 km, and A=1.0 cm, 



0.67 x 10 25 7 cm . , 
Pns = M ns v ns c — — T 9 gm— , (16) 

with Tg = 1( ^ K ■ Taking the mass of the neutron 
star to equal the mass of our sun, M ns = 2 x 10 33 
gm, we obtain for the velocity of the neutron star 

-«5X10-(^_1^^,(1 7 ) 



For example, for T=10 Mev =1.16 x 10 11 kK, 

v na ~ 47.3— to 2, 370— , (18) 
s s 

which means that sterile neutrino emission could ac- 
count for the large pulsar kick with the parameters 
extracted from Refs [1, Q. If we use the physical 
parameters that give Eq. (fl8)) for electron neutrinos, 
we obtain a pulsar velocity of v ns — 95 km/s, which 
is consistent with previous predictions by several au- 
thors. 

It should be noted that the study of the Mini- 
BooNE and LSND results are in progress, and the 
mixing angles that result could be much different 
from those which we have used, changing our re- 
sults. Preliminary data from the MiniBooNE/Minos 
experiment [25[ is consistent with the MiniBooNE 
results[7J. We also once more point out that Ref [10( 
questions the accuracy of the peramaters extracted 
by Ref ||. 
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III. CONCLUSIONS 

Because of the strong magnetic fields in protoneu- 
tron stars and the associated neutrinosphere, the 
electrons produced in the URCA processes that 
dominate neutrino production in the first fO sec- 
onds have a sizable probability, P(0), to be in the 
lowest (n=0) Landau level. This leads to asymmet- 
ric neutrino momentum. With the mixing angles 
found in Refs[l, Q, we find that the sterile neutri- 
nos produced during this period for high luminos- 
ity pulsars can give the pulsars velocities of greater 
than 1000 km/s, as observed, similar to predictions 
based on sterile neutrinos as dark matter M. Al- 
though the parameters found in the analyis [H, Q of 
the MiniBooNE Q and LSND data is still in question, 
current analysis of the FermiLab experiment [25j 
seems to be in agreement with MiniBooNE 0]. 

There is a strong correlation of the pulsar velocity 
with temperature, T. Since it is difficult to determine 
T accurately, it is difficult for us to predict the ve- 



locity of a pulsar whose kick arises from sterile neu- 
trino emission. On the other hand, if the pulsar kick 
arises from the asymmetric emission of active neutri- 
nos produced by the modified URCA processes after 
10 seconds, also proportional to P(O)0], then T can 
be determined by an accurate measurement of the 
neutrinos from the supernova. Therefore, in future 
years, with much more accurate neutrino detectors, 
one could predict the velocity of the resulting pulsar. 
Unfortunately, the energy of emitted sterile neutri- 
nos cannot be measured. From our results in the 
present paper and those in RefQ, high luminosity 
pulsars receive a large kick both from sterile neutri- 
nos in the first ten seconds and standard neutrinos 
in the second ten seconds. 
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